function [m] = min_distance(est_params,glob,param)

    %{
    This function simply defines the minimum distance 
    criterion. (forces it to evaulate solve_simulate
    only once per iteration)
    %}

    dist = solve_simulate(est_params,glob,param,-1,0) - glob.moments_data;
    m = dist'*inv(glob.W)*dist;


end